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Abstract 

History, status and plans of the search for critical behavior of strongly interacting matter created 
in nucleus-nucleus collisions at the CERN Super Proton Synchrotron is reviewed. In particular, it is 
expected that the search should answer the question whether the critical point of strongly interacting 
matter exists and, if it does, where it is located. 

First, the search strategies are presented and a short introduction is given to expected fluctuation 
signals and to the quantities used by experiments to detect them. The most important background 
effects are also discussed. 

Second, relevant experimental results are summarized and discussed. It is intriguing that both 
the fluctuations of quantities integrated over the full experimental acceptance (event multiplicity and 
transverse momentum) as well as the bin size dependence of the second factorial moment of pion and 
proton multiplicities in medium-sized Si-|-Si collisions at ISSvl GeV/c suggest critical behaviour of the 
created matter. 

These results provide strong motivation for the ongoing systematic scan of the phase diagram by 
the NA61/SHINE experiment at the SPS and the continuing search at the Brookhaven Relativistic 
Hadron Collider. 

PACS numbers: 12.40.-y, 12.40.Ee 
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I. INTRODUCTION 


The structure of the phase diagram of strongly interacting matter is one of the most impor¬ 
tant topics in nuclear and particle physics. We know that at low densities strongly interacting 
particles are hadrons and thus the matter is in form of a hadron gas or liquid. Since the discov¬ 
ery of sub-hadronic particles, quarks and gluons, it was speculated that at high temperature 
and/or pressure densely packed hadrons will ’’dissolve” into a new phase of quasi-free quarks 
and gluons, the quark-gluon-plasma (QGP) [1]. 

Many years of intense experimental and theoretical studies of high energy nucleus-nucleus 
(A-l-A) collisions led to the conclusion that the quark-gluon plasma exists in nature. This 
conclusion is based on a wealth of systematic data on A-l-A collisions at very high energies 
from the CERN Large Hadron Collider (LHC) [2] and the BNL Relativistic Heavy Ion Collider 
(RHIC) (see e.g. Ref. [3]), and, very importantly the observation of the transition between 
hadronic matter and quark-gluon plasma at the CERN Super Proton Synchrotron energies [1] 
(for recent review see Refs. [SI [6]). 

Thus the current key question in the study of the phase diagram of strongly interacting 
matter is the structure of the phase transition region between the hadron gas and the quark- 
gluon plasma. 

The most popular suggestion is shown schematically in Fig. The transition at low tem¬ 
perature T and high baryochemical potential /ib is believed to be of the first order and happen 
along a line which ends with decreasing /xb in a critical point (of the second order) and then 
turns into a crossover region. The non-trivial structure of the phase transition region was first 
suggested by Asakawa, Yazaki [7] and Barducci, Casalbuoni, De Curtis, Catto, Pettini [S]. Ex¬ 
perimental studies of the features of the phase diagram were strongly motivated by predictions 
of measurable effects. The pioneering work of Wosiek [9] pointed out that intermittent [10] 
behaviour is naturally expected at a phase transition of the second order. Soon after the con¬ 
jecture was further developed by Satz m, Antoniou et ah na and Bialas, Hwa [13]. This 
initiated experimental studies of the structure of the phase transition region via studies of 
particle multiplicity fluctuations using scaled factorial moments. Later additional measures 
of fluctuations were also proposed as probes of critical behaviour [sms]. The results of the 
KLM [16] and NA49 [IT] experiments suggest that effects related to the critical point may have 
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been observed in collisions of medinm size nuclei at the top SPS energy. This motivated the 
NA61/SHINE experiment to perform a systematic scan in collision energy and system size. The 
new measurements should answer the general question about the nature of the transition region 
and, in particular, the question; does the critical point of strongly interacting matter 
exist in nature and, if it does, where is it located? The most recent experimental and 
theoretical status of the exploration of the phase diagram is reviewed at the regular workshops 
on the Critical Point and Onset of Deconhnement na. 



Figure 1: Sketch of the phase diagram of strongly interacting matter. Italic labels added to the 
sketch show regions probed by the early stage of heavy ion collisions studied in current (in red) and 
future (in green) experimental programmes. 

The structure of the transition region is explored experimentally by studying the hnal states 
produced in nucleus-nucleus collisions. By changing collision energy and size of colliding nuclei 
one changes temperature T and chemical potentials /is of matter at the freeze-out stage [T^ . 
In particular, by increasing collision energy more and more pions per colliding baryon are 
produced which is the main cause for the decrease of the baryon chemical potential with collision 
energy. With increasing size of the colliding nuclei the volume of created matter increases 
and consequently the role of hadron-hadron interactions at the late stage of matter evolution 
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Figure 2: Left: Compilation of chemical freeze-out points of central Pb-|-Pb (Au-|-Au) collisions [20]. 
Right: Chemical freeze-out points of nucleus-nucleus collisions studied with the NA49 |19) and 
NA61/SHINE |3T] programmes at the CERN SPS. 

becomes more important. This leads to a decrease of the freeze-out temperature with increasing 
size of the colliding nuclei. Thus scanning in collision energy and system size one hopes to be 
able to move the freeze-out close to the transition region. This is illustrated in Fig. The 
left plot shows the collision energy dependence of the freeze-out parameters in central Pb-|-Pb 
collisions [2TH23], whereas the right one presents their energy and system size dependence at 
the CERN SPS [I9]. 

The experimental search for the critical point by investigating nuclear collisions is promising 
only at energies higher than the energy of the onset of deconhnement, which experimentally 
was located to be at the low SPS energies [211125]. This is because the energy density at the 
early stage of the collision, which is required for the onset of deconhnement is higher than the 
energy density at freeze-out, which is relevant for the search for the critical point. 

A characteristic feature of a second order phase transition (the critical point or line) is 
the divergence of the correlation length. The system becomes scale invariant. This leads to 
large huctuations in particle multiplicity. Moreover these huctuations have specihc characteris¬ 
tics mm- Also other properties of the system should be sensitive to the vicinity of the critical 
point [15]. Thus when scanning the phase diagram a region of increased fluctuations may signal 
the critical point or the critical line. This is illustrated schematically in Fig. [^ which presents 
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an updated version of the plot shown hrst in Ref. [^B] . 



Figure 3: Sketch of the ’’critical hill” expected in the search for the critical point in the two dimen¬ 
sional plane (system size) - (collision energy). At the hill the characteristic fluctuation signals of the 
critical point are maximal, see Sec. ra for details. 


The study of fluctuations and correlations is significantly more difficult than the study of 
single particle spectra and mean multiplicities. In general, results on fluctuations are sensitive 
to conservation laws, resonance decays and many of them also to the unavoidable volume 
fluctuations of colliding nuclear matter. Moreover, they cannot be corrected for a limited 
experimental acceptance. 

This review is organized as follows. In Sec. [IT| experimental strategies, as well as techniques 
and problems for the search for the critical point are briefly presented. Search results from 


experiments at the CERN SPS, in particular NA49 and NA61/SHINE, are reviewed in Sec. Ill 


Conclusions and an outlook in Sec. IV close the paper. 
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II. SEARCH STRATEGIES, TECHNIQUES AND PROBLEMS 


This section reviews basic ideas and tools relevant for the experimental search for the critical 
behaviour of strongly interacting matter at the CERN SPS. The most important background 
effects are listed and examples are discussed. 


A. Onset of deconfinement versus critical point 


Onset of Deconfinement: 



Critical Point: 




Figure 4: Sketch of the search strategy for the onset of deconfinement (left) and the critical point 
(right). The location in the phase diagram of the matter at the early stage is indicated by the red dot 
and at freeze-out by the blue square. The solid (dotted) lines show the first order (rapid crossover) 
boundary and the critical point is the grey dot at the end of the first order transition line. 
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Here relations between the onset of deconfinement, the critical point of strongly interacting 
matter and the possibilities of their experimental study in relativistic nucleus-nucleus collisions 
are discussed. The two sketches presented in Fig. [^should help to understand the basic ideas. 

The onset of deconhnement refers to the beginning of the creation of a deconhned state of 
strongly interacting matter (ultimately a quark-gluon plasma) at the early stage of nucleus- 
nucleus collisions when increasing the collision energy. With increasing collision energy the 
energy density of matter created at the early stage of A-l-A collisions increases^. Thus, if there 
are two phases^ of matter separated by the transition region (solid and dotted lines) as indicated 
in Fig. 1^ left the early stage (the red point) has to cross the transition region. Therefore, the 
existence of the onset of deconhnement is the most straightforward consequence of the existence 
of two phases of strongly interacting matter, i.e. conhned matter and QGP. The experimental 
observation of the onset of deconhnement required a one dimensional scan in collision energy 
with heavy ions as performed by NA49 [21]. Signals of the onset of deconhnement relate to the 
diherence in properties of conhned matter and QGP. They are weakly sensitive to the structure 
of the transition region. 

Discovery of the onset of deconhnement implies the existence of QGP and of a transi¬ 
tion region between conhned and QGP phases. Recent experimental results [21| indicate that 
the transition region (the mixed phase or rapid cross-over) ranges from ^/sNN ~ 8 GeV to 
^/sNN ~ 12 GeV, where ^v/snn denotes collision energy per nucleon pair in the centre of mass 
system. 

Numerous possibilities concerning the structure of the transition region are under discussion 
(see e.g.. Ref. my The most popular one [28l [29|, sketched in Fig. claims that a 1'^* order 
phase transition (thick gray line) separates both phases in the high baryonic chemical potential 
domain. In the low baryonic chemical potential domain a rapid crossover is expected (dotted 
line). The end point of the order phase transition line (grey dot in Fig. is the critical 

^ The correlation between collision energy and the early stage energy density is expected to be strong for 
central collisions of large nuclei and weak for collisions of low mass nuclei. In the latter case energy and/or 
multiplicity of produced particles may help select collisions with a similar energy density. 

^ The discussed two phase diagram is the simplest one which allows to introduce the concepts of the onset 
of deconhnement and the transition region. There are numerous suggestions of phase diagrams with a much 
richer structure (see e.g., Ref. m)- 






point (of the second order). 

The characteristic signatnres of the critical point can be observed if the freeze-out point 
(blue square in Fig. right) is located close to the critical point. The analysis of the existing 
experimental data [19] indicates that the location of the freeze-out point in the phase diagram 
depends on the collision energy and the mass of the colliding nuclei. This dependence is 
schematically indicated in Fig. right and quantihed in Fig. right. Thus the experimental 
search for the critical point requires a two-dimensional scan in collision energy and size of the 
colliding nuclei. The NA61/SHINE experiment [26l[3T] at the CERN SPS started this scan in 
2009 and completion is expected within the coming few years. Note, that a two dimensional 
scan is actually required for any study of the structure of the transition region, independent of 
the hypothesis tested. 

The transition region can be studied experimentally in nucleus-nucleus collisions only at T, 
fiB values which correspond to collision energies higher than the energy of the onset of decon- 
hnement. This important conclusion is easy to understand when looking at Fig|^ Signals of the 
transition region can be observed provided the freeze-out point is close to it (see Fig. |^ng'/if). 
Furthermore, the energy density at the early stage of the collision is, of course, higher than the 
energy density at freeze-out. Thus, the condition that the freeze-out point is near the transition 
region implies that the early stage of the system is above (or on) it. This in turn means that 
the optimal energy range for the search for the critical point (or, in general, for the study of 
properties of the transition region) lies above the energy of the onset of deconfinement (see 
Fig. 1^ left). This general condition limits the search for the critical point to the collision energy 
range Elab > 30A GeV {^/sNN ~ 8 GeV). 

B. Scaled factorial moments 

In the grand canonical ensemble the correlation length ^ diverges at the critical point (or 
second order phase transition line) and the system becomes scale invariant [H E]. This leads 
to large multiplicity fluctuations with special properties. They can be conveniently exposed 
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using scaled factorial moments Fr{6) [10] of rank (order) r: 


FAS) = ^ -. ( 1 ) 

i=l 

where 6 is the size of the subdivision intervals of the momentum phase space region A and 
M = A/5 is the number of intervals. Ni refers to particle multiplicity in the interval i and (...) 
indicates averaging over the analyzed collisions. 

For a non-interacting (ideal) gas of Boltzmann particles in the grand canonical ensemble 
(IB-GCE) one gets Fr{6) = 1 for all values of r and 6 provided the mean particle multiplicity is 
proportional to 6. The latter condition is trivially obeyed for a subdivision in configuration space 
where the particle density is uniform throughout the gas volume. For the case of subdivision 
in momentum space the subdivision should be performed using so-called cumulative kinematic 
variables [32] in which the particle density is uniform. 

At the second order phase transition the matter properties strongly deviate from the ideal 
gas. The system is a simple fractal and F)(5) possess a power law dependence on S: 

Fri6) = FriA) ■ {A/5f^ . (2) 

Moreover the exponent (intermittency index) satishes the relation: 

(pr = {r -1) ■ dr , (3) 

with the anomalous fractal dimension dr being independent of r [TB]. These results are valid 
when cumulative variables [32] are used to define the intervals 6. The properties represented 
by Eqs. [^ and are called in this paper the critical behaviour of the scaled factorial moments. 

An experimental search for the properties[^and[^in high energy collisions requires signihcant 
additional input. In particular, one has to decide on: 

(i) dimension, size and location of the momentum phase space region A, 

(ii) selection of collisions used in the analysis, 

(iii) selection of particles used in the analysis. 
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Concerning (i) it was shown by Bialas, Seixas [33] (see also Ochs, Wosiek [331 I3S]) that 
nnbiased results on the critical behaviour of scaled factorial moments can be obtained only by 
performing the analysis in variables and dimensions in which the singular behaviour appears. 
Any projection procedure is likely to to remove, at least partly, the critical fluctuation signal. 

Concerning (iii) QCD-inspired considerations [2HI ES] suggest that the order parameter of 
the phase transition is the chiral condensate {qq) {q is the quark held). The quantum state 
carrying the quantum numbers as well as the critical properties of the chiral condensate is the 
isoscalar cr-held. Assuming that this state can be formed in high energy collisions there are two 
possibilities for the observation of the properties and 

(i) Directly from its decay products 1371 . The condensate will decay into tt^tt pairs with 
invariant mass just above twice the pion mass. Detection of the expected huctuations 
requires reconstruction of the pion pairs. One expects here d = 02 = 2/3 [37] . 

(ii) Through measuring the huctuations of the proton number. The net-baryon density mixes 
with the chiral condensate transferring the critical huctuations to the net-baryon den¬ 
sity [28l [381 - 133] . which is an equivalent order parameter of the phase transition. The re¬ 
sulting huctuations are predicted to be present also in the net-proton number as well as in 
the proton and anti-proton numbers separately [13]. One expects here d = 02 = 5/6 |45j . 


C. Central moments 

The inhnite correlation length at a second order phase transition is expected to lead to 
divergence of the second central moment of the multiplicity distribution. This was recently 
illustrated by Vovchenko, Anchinshkin, Gorenstein and Poberezhnyuk [361I3Z] via analytical 
calculations performed in the GCE for a gas obeying the van der Waals equation of state. As 
seen in Fig. |^the scaled variance oo = {{N‘^) — {N)‘^)/{N) increases when the critical point is 
approached and it diverges at the critical point. 

Moments of other quantities are also expected to increase in the vicinity of the critical 
point. In particular, these are event transverse momentum (vector magnitudes are summed) 
and net charges of produced particles [13]. Of course the latter have to be studied in a phase 
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Figure 5: Left: Scaled variance of the multiplicity distribution calculated within the grand canonical 
formulation of the van der Waals model (massless particles obeying Boltzmann statistics) as a function 
of temperature T and particle density p scaled by the corresponding values at the critical point (from 
Ref. |36] ). Right: Scaled variance of the multiplicity distribution calculated within the grand canonical 
formulation of the van der Waals model (nucleons obeying Fermi statistics) as a function of temperature 
T and baryon chemical potential p, (from Ref. [47jb 


space region which covers only a small fraction of all particles in order to avoid suppression of 
fluctuations due to conservation laws (see below for details). 

It was pointed out by Stephanov [18] that second moments of the multiplicity distribution 
increase in proportion to the square of the correlation length while moments of higher order 
are proportional to even higher powers of the correlation length. Moreover, ratios of certain 
higher order cumulants are expected to be independent of the correlation length and therefore 
allow a consistency test Based on QCD model calculations the same paper also presents 
quantitative estimates of the size of the enhancement of multiplicity fluctuations caused by a 
critical point. 


D. Strongly intensive quantities 

Since event-to-event volume fluctuations cannot be eliminated in experimental studies of 
nucleus-nucleus collisions, it is important to minimise their effect by dehning suitable fluc¬ 
tuation measures. It appears that using second and hrst moments of the distibution of two 
extensive quantities (their hrst moments are proportional to volume) one can construct huctu- 
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ation measures which are, within a statistical model of the ideal Boltzmann gas in the grand 
canonical ensemble formulation ( SM(IB-GCE) ) |19], independent of its fluctuations. 

As the simplest example let us consider multiplicities of two different types of hadrons, A 
and B. Their mean multiplicities are proportional to the system volume: 

{A) r^V , {B) r^V . (4) 

Obviously the ratio of mean multiplicities is independent of the volume V. Moreover, the ratio 
(A) / {B) is independent of PiV), where PiV) is the probability (density) distribution of V for 
the considered set of collisions. Quantities which have the latter property are called strongly 
intensive quantities [19]. Such quantities are recommended to be used in experimental studies 
of the system size dependence of fluctuations in A+A collisions as they eliminate the influence 
of usually poorly known distributions of the system volume. 

More generally, A and B can be any extensive event quantities such as the sum of transverse 
momenta, the net charge or the multiplicity of a particular type of particle. The scaled variances 
of A and B and the mixed second moment {AB) calculated within the SM(IB-GCE) [35] read: 


{AB) 


u[A]=u*[A] + {A)/{V)-u[V] , 

(5) 

u:[B]=u*[B] + {B)/(y)-u:[V] , 

(6) 

{ABY{V) + {A){B){V)P{{V^)-{V)) , 

(7) 


where quantities denoted by * are quantities calculated for a fixed value of the system volume. 

From Eqs. [sjj^ follows jlHl EQ] that suitably constructed functions of the second moments, 
namely 


A[A,B] 


{B)u;[A] - {A)u[B] 


( 8 ) 


and 


E|AB] = 7=r + (.4 )i.,[B] - 2{{AB) - {A){B)) 

Os L 


(9) 


are independent of P{V) in the SM(IB-GGE). Here the normalisation factors Ca and Gs are 
required to be proportional to first moments of any extensive quantities. In Ref. |SI a specific 
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choice of the Ca and Cy, normalisation factors was proposed which makes the qnantities A [A, B\ 
and E[A,i?] dimensionless and leads to A[A,B] = Ti[A,B] = 1 in the independent particle 
model (IPM). This normalisation is called here the IPM normalisation and unless otherwise 
stated the IPM normalisation is used. 

Thus A [A, B] and S[A, B] are strongly intensive quantities which measure fluctuations of A 
and B, i.e. they are sensitive to second moments of the distributions of the quantities A and B. 
Results on A [A, B] and E[A, B] are referred to as results on A — B fluctuations, e.g. transverse 
momentum - multiplicity fluctuations. 

For the case of multiplicity A - multiplicity B fluctuations the expressions for A[A,B] and 
Tj[A,B] with the IPM normalisation read: 

A[A, B] = {{B)u:[A] - {A)u:[B]) / {{B) - {A)) (10) 

and 

T.[A,B] = ({BVM + - 2((AB) - {A){B))) / ((B) + {/1» . (11) 

The S quantity is a reincarnation of the popular $ measure of fluctuations isj. The original 
dehnition of $ is: 



where x is a particle property and 

N 

Z = ^{xi-x) ( 13 ) 

i=l 

with the sum running over the N particles in the event. $ was shown to be related to E as 
follows [i9] : 

<I>a; = aJxu:[x\ 

with X = = {N)uj[x]. 

By construction strongly intensive measures of fluctuations are always a functional of two 
extensive quantities. This, in general, hampers a straight-forward interpretation of experimental 
results. However, under certain conditions the A quantity can be used to obtain the scaled 
variance of the extensive quantity A separately. 

Let A be the extensive quantity, e.g. selected for its sensitivity to the critical fluctuations. 
Then choose a quantity B such which is proportional to the system volume B ^ V. Then it 


VS[A,A]-1 


( 14 ) 
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is easy to show that the strongly intensive measures Ab[v4, B] and B] (equal to A [A, B] 

and i?] with the normalisation Ca = {B) ~ V) obeys the relation: 

Ab[A,B] = J:b[A,B]=u*[A] . (15) 

In the derivation of Eq. one assumes the validity of Eq. which needs to be investigated 
case-by-case. Thus Ab[A,B] is approximately equal to the scaled variance Ci;*[A] of A for a 
hxed system volume (see Eq. : 

Ab[A,Bv]^uj*[A] , (16) 

Suggestions of practical choices of A and B are: 

(i) A - multiplicity of hadrons which are sensitive to the critical behaviour, e.g. protons, 
sigma-mesons, pions in the central rapidity window and 

(ii) B - net electric charge in full phase space or large acceptance excluding fragmentation 
regions of projectile and target nuclei; it is equal to the number of participant protons 
and thus approximately proportional to the volume of matter involved in the collision or 
- the number of projectile participants calculated as the difference between the number of 
nucleons in the beam nuclei and the number of projectile spectators measured by a ’’zero 
degree” calorimeter, e.g., the Projectile Spectator Detector of NA61/SHINE. 

Recently, strongly intensive measures which involve higher than second moments were pro¬ 
posed |S2]- The next important step would be to reformulate the critical properties of scaled 
factorial moments Eqs. |^and[^in terms of strongly intensive quantities. 


E. Critical hill 


As discussed in Sec. |II C| moments of multiplicity distribution diverge when temperature and 
density approach their critical values (see Fig. [^. At this same point in the phase diagram the 
scaled factorial moments should obey the critical properties described by Eqs. and 

As previously argued the freeze-out temperature increases with increasing collision energy 
and decreasing size of the colliding nuclei. The maximum temperature, probably the closest 
to the phase transition, is observed in p-|-p interactions. However, the small volume and short 


15 



life-time of the created matter together with the conservation laws (see the next section), do not 
allow the divergence of the correlation length [53l IS]. Thus neither a divergence of fluctuation 
measures, e.g. of the scaled variance, nor the appearance of critical behaviour of scaled factorial 
moments are expected. 

On the other hand, critical fluctuations developing in the system produced by collisions of 
large nuclei may be erased by re-scattering processes between hadronisation and kinetic freeze- 
out. Therefore one can argue that the maximal signals of the critical point may be observed for 
collisions of medium mass nuclei. In this case the volume and life-time of the created matter 
are large enough to allow for the critical behaviour to appear, and the temperature is close 
enough to the critical temperature to make the critical behaviour visible. Thus observation of 
the ’’critical hill”, as sketched in Fig. would provide convincing evidence for the existence 
and location of the critical point. 


F. Background fluctuations 

The critical hill, if observed, will rise above a background caused by many different sources 
of fluctuations. These are in particular: 


(i) volume fluctuations discussed in Sec. IID 


fii) conservation laws discussed in Sec. IIG 


(hi) formation and decay of resonances, 

(iv) quantum statistical effects (Bose-Einstein and Fermi-Dirac statistics). 

Their impact on fluctuations is discussed in the recent review [55] where also references to 
original papers are given. 

In general the background fluctuations are not expected to lead to non-monotonic depen¬ 
dence of fluctuations on collision energy and system size. 


The effect of volume fluctuations is addressed in Sec. IID As an additional example, the 
effect of conservation laws is discussed below in more detail. 
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G. Conservation laws 


Predictions of statistical models concerning the volume dependence change qualitatively 
when material and/or motional conservation laws are introduced, i.e. instead of the grand 
canonical ensemble, the canonical (CE) or micro-canonical (MCE) ensembles are used. The 
effect of conservation laws has been extensively studied for mean multiplicities since 1980 (see 
e.g. Refs. [SnHSH]) and for second moments of multiplicity distributions since 2004 (see e.g. 
Refs. [HHIEO]). An example is discussed below for illustration. 

Figure]^ taken from Ref. [59] presents the results of calculations performed within the sim¬ 
plest model which allows to study the effect of material conservation laws on mean multiplicity 
and scaled variance of the multiplicity distribution. In this model an ideal gas of classical 
positively and negatively charged particles is assumed. The ratio of the mean multiplicities 
calculated within the SM(IB-CE) and the SM(IB-GCE) is plotted in Fig. |^(/e/f) as a function 
of the mean multiplicity 2 ; from the SM(IB-GCE), the latter being proportional to the system 
volume. The ratio approaches one with increasing volume. Thus for sufficiently large systems 
mean multiplicities obtained within the SM(IB-GGE) can be used instead of mean multiplicities 
from the SM(IB-GE) and the SM(IB-MGE) [60]. This is however not the case for the scaled 
variance as illustrated in Fig. [right). The results for the SM(IB-GE) and the SM(IB-GGE) 
approach each other when the volume decreases to zero. Of course, the scaled variance in 
the SM(IB-GGE) is one independent of volume. Different behaviour is observed for the scaled 
variance in the SM(IB-GE), where it decreases with increasing volume and for a sufficiently 
large volume approaches 0.5. 

The above example clearly illustrates the importance of material and motional conservation 
laws for hadron production in nucleus-nucleus collisions at high energies. In particular, the 
conservation laws may strongly affect ffuctuations even for large systems, leading to signihcant 
deviations from the simplest reference model, i.e. independent particle production. 

III. SEARCH RESULTS FROM EXPERIMENTS AT THE CERN SPS 

This section reviews the status of the experimental search for evidence of a second order 
phase transition and/or the critical point at the GERN SPS based on published results and 
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Figure 6: The ratio of mean multiplicities {N) = (A^+) or {N-) calculated within the SM(IB-CE) 
and the SM(IB-GCE) are plotted in the left panel as a function of mean multiplicity z for the SM(IB- 
GCE), with z being proportional to the system volume. The corresponding scaled variances w'*', u}~ are 
shown in the right panel. The ideal gas model of classical positively and negatively charged particles 
was used for calculations. The system net charge is assumed to be zero and thus (A^+) = {N-) and 
= {uj~). The plots are taken from Ref. [59] . 

preliminary data presented at conferences. 


A. Pioneering analyses 

The search for the critical behaviour of strongly interacting matter at the CERN SPS started 
in 1990 from the pioneering paper of Bialas and Hwa [13] . In this work the authors compiled 
results from intermittency analysis performed by the EMC [61], NA22 [6^ and KLM [T6] 
experiments at the SPS. In these experiments intermittency indices (pr were determined from 
power-law hts to the dependence on bin-size in (pseudo-)rapidity of the scaled factorial moments 
of successive rank r. Figure [^shows that the anomalous dimensions, defined as dr = (pr /~ 1); 
increase strongly with rank r for the smaller reaction systems. Interestingly the anomalous 
dimension stays constant for the heaviest system, for S-|-(Ag Br) collisions at 200A GeV/c. 
Bialas and Hwa interpreted this behaviour as an indication of a second order phase transition. 
A quantitative prediction of dr was not provided. 
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Figure 7: Anomalous dimensions dr = <^r/{T — 1) of pseudo-rapidity spectra of hadrons produced in 
/i-|-p, h-|-p, p-|-AgBr, 0-|-AgBr and S-|-AgBr collisions at ^ 20 GeV [13]. 

An extended analysis of the emulsion data of the KLM collaboration was published in 
Ref. [nSj- Corrections were applied for the non-uniform rapidity distribution and the inter- 
mittency indices were determined both for 1-dimensional (pseudo-rapidity) and 2-dimensional 
(pseudo-rapidity and azimuthal angle) subdivisions of phase space. The resulting anomalous 
dimensions for S-|-(Ag Br) collisions at 200^4 GeV/c are plotted in Fig. One observes that 
the values of dr are consistent with being independent of r for both analyses and conhrm the 
earlier results. However, the values of dr are roughly 5 times larger in the 2-dimensional analy¬ 
sis. A strong reduction of the measured power (pr with decreasing dimensionality of the analysis 
was explained by Bialas and Seixas [33] as due to averaging of fluctuations via the projection 
procedure. Thus the factorial moment analysis in three dimensions seems to be mandatory in 
future searches for the critical behaviour. 

Motivated by these results the WA80 [M] and NA35 [65] experiments at the SPS revisited 
intermittency analysis in nucleus-nucleus collisions at 200A GeV/c. WA80 did not have mo¬ 
mentum measurement and inferior 2-track and angular resolution compared to the emulsion 
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Figure 8: Left: Scaled factorial moments of rank i from two-dimensional analysis in pseudo-rapidity 
and azimuthal angle as a function of subdivision size for the 19% most central S-|-(Ag Br) collisions 
at 200A GeV/c [63]. Right: The corresponding anomalous dimensions dr versus order r are shown by 
squares, whereas the circles show results for the one-dimensional analysis in pseudo-rapidity. 

experiment. They concluded that they observed no signihcant intermittency effect in S-l-S and 
S-I-Au collisions when taking into account statistical and systematic uncertainties. 

The NA35 streamer chamber experiment performed momentum measurements in central 
p-|-Au, 0-|-Au, S-l-S and S-|-Au collisions which were subjected to a fully differential 3- 
dimensional (rapidity, transverse momentum, azimuthal angle) intermittency analysis. Al¬ 
though the factorial moments were found to increase with the number of subdivisions of phase 
space this rise was not well described by a power law. Instead, a conventional model supple¬ 
mented by Bose-Einstein correlations provided a satisfactory description. 

In summary, the intermittency analyses of charged particle production in oxygen and sulphur 
induced reactions did not lead to conclusive results on the existence of a second-order phase 
transition in these reactions. More recent theoretical investigations suggest that when the 
hadronization of a QGP occurs near the critical point the hadronization of the chiral condensate 
will lead to intermittency in the production of protons and low-mass tt+tt" pairs with known 
intermittency index. A search for such effects is in progress in the NA49 experiment and will 
be discussed below. 
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Figure 9: Reactions and energies of the ongoing scan of the phase diagram by NA61 {left) (filled green 
squares denote completion of data taking) and systems previously studied by NA49 {right). 

B. Systematic fluctuation studies of NA49 and NA61/SHINE 


Search for evidence of critical behaviour of strongly interacting matter was restarted at the 
beginning of 2000 by the NA49 Collaboration [66]. Fluctuations were analyzed in a large number 
of previously recorded data sets (see Fig. right). A tantalizing increase of multiplicity and 
transverse momentum fluctuations was found in collisions of medium size nuclei at 158A GeV/c. 
This motivated the ongoing measurements of NA61/SHINE [67], the successor of the NA49 
experiment. For the hrst time a systematic two-dimensional scan in system size and collision 
energy is being performed (see Fig. |^Ze/f). A search for the critical point is also in progress 
at RHIC by the STAR collaboration [3] . Relevant results will be mentioned at the appropriate 
places. The results from the SPS obtained by NA49 and NA61/SHINE are reviewed below. 

When comparing results from the SPS to those from the RHIC-BES program one should 
keep in mind that the event selection and acceptance of the experiments are signihcantly dif¬ 
ferent. NA49 [HH] and NA61/SHINE jUl] are hxed-target spectrometers and cover mainly the 
forward region in the center-of-mass system, particularly when particle identihcation is re¬ 


quired (see Eig. 10). An advantage of the hxed-target geometry is that it allows to characterise 
the centrality of the collisions by measuring the energy of the spectators from the beam nu¬ 
cleus independently from the measurements performed on the produced particles. On the other 
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Figure 10: Acceptance of the NA49 and NA61/SHINE experiments for pions (a), kaons (b) and protons 
(c) at the SPS at 30A GeV (geometrical acceptance in grey, acceptance for identification in color), as 
well as that of the STAR detector at RHIC (lines) at the equivalent energy ^sjviv = 7.7 GeV. 


hand, STAR at RHIC is a collider experiment with practically energy independent rapidity ac¬ 


ceptance \y\ < 0.7, but without the low transverse momentum region (see curves in Fig. 10). 
The track density in the detector increases only moderately with collision energy. However, 
the projectile spectator regions are not accessible to measurement and the collision centrality 
selection has to be based on the multiplicity of produced particles. 
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Figure 11: Second scaled factorial moments F 2 {M) of the proton number in transverse momentum 
space at mid-rapidity (—0.75 < y < 0.75) for the most central collisions of (a) “G”-)-C (12%), (b) 
“Si”-|-Si (12%), and (c) Pb-|-Pb (10%) at ^smn = 17.3 GeV. The circles (crosses) represent F 2 {M) of 
the data (mixed events) respectively. 
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1. Scaled factorial moments 


The NA49 experiment at the CERN SPS searched for an intermittency signal in the pro¬ 
duction of proton EH and low-mass vr^'^Tr pairs Ea in the most central collisions (12%, 12%, 
10%) of “C”, “Si” and Pb nuclei on C (2.4% interaction length). Si (4.4%) and Pb (1%) tar¬ 
gets, respectively, at beam energy of IhSA GeV (yGiviv = 17.3 GeV). The analysis looked in 
transverse momentum space for a power law behaviour of the second scaled factorial moments 
(SSFMs, F 2 {M)) defined in Eq. [^with r = 2. 

4.1 Intermittency in proton production at mid-rapidity 

Protons were identified with a purity of above 80% based on the ionization energy loss of 
the tracks in the TPG detectors. Poorly measured and fake tracks were carefully removed 
by suitable selection criteria since they can produce a spurious intermittency signal. Since 
critical fluctuations originating from the GP are predicted to be strongest in a region around 
mid-rapidity in the cms system, protons were selected in the rapidity range —0.75 < y < 0.75. 

The strong background from mis-identified and non-critical protons was estimated and sub¬ 
tracted using mixed events, which by construction do not contain critical fluctuations. The 


resulting SSFMs are plotted in Fig. 11 Evidently the mixed event background is consistent 
with the data for “G”-l-G and Pb-I-Pb collisions. On the other hand, in “Si”-1-Si reactions the 
values of F 2 {M) increase with while those of the mixed events remain nearly constant, 
shows the background subtracted SSFMs F!f'\M) 


Figure 


12 


AF^^\m) = F:f\M) - Fr\M) , 




{m)i 


(17) 


where F!f\M) and F!f^\M) correspond to the data and mixed-event background, respec¬ 
tively. While the results for “G”-l-G and Pb-|-Pb collisions scatter around zero, the values for 
“Si”-|-Si reactions rise with M^. A power-law fit in the region > 6000 gave the result 
(f )2 = 0.96lo;25(stot) ± 0.1Q{syst) with y^/dof 0.09 - 0.51. The re-sampling method [73] was 
used to calculate the errors in order to take account of the strong correlation between successive 
(in M) values of the SSFMs. Thus no intermittency signal is present for proton production in 
“G”-|-G and Pb-I-Pb collisions whereas power-law fluctuations are observed in “Si”-|-Si reactions 
with an exponent consistent with the GP prediction. 

4-2 Intermittency in low-mass 7 r~^n~ pair production 
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Figure 12: The combinatorial background subtracted moments /S.F!f\M) corresponding to the mo¬ 


ments of Fig. 11 in bins of transverse momentum for the most central collisions of (a) “C” +C (centrality 
12 %), (b) “Si”+Si (centrality 12 %) and (c) Pb+Pb (centrality 10 %) at = 17.3 GeV. The line 

in the middle plot shows the result of a power-law fit for > 6000 with exponent 0.96. 


Results of a search for critical fluctuations in the chiral condensate via a similar intermittency 
study of low-mass tt+tt" pairs was published by the NA49 collaboration in Ref. |7Z|. The chiral 
condensate is believed to decay into tt+tt” pairs near the mass threshold when deconfined matter 
hadronises. As in the case for protons these fluctuations may be detectable by studying SSFMs 
of the TT+TT" pair number, provided the combinatorial background can be sufficiently reduced. 

Pions were required to have laboratory momenta exceeding 3 GeV/c and identification was 
based on the ionization energy loss of the tracks in the TPC detectors. For further analysis 
low-mass tt+tt" pairs satisfying 


2mTr + ei < m^r+TT- < 2777-,^ + 62 


(18) 


were considered, where is the pair invariant mass, ei = 5 MeV was chosen to remove the 

enhancement of pairs from Coulomb attraction and 62 = 34,24,1 MeV for “C”-|-C, “Si”-|-Si, 
and Pb-|-Pb respectively was optimised to reduce the combinatorial background. Finally the 
remaining background was estimated by pairs from mixed events which were made to satisfy 
the same criteria. 

The SSFMs of the tt+tt” pair multiplicity distribution versus the number of subdivisions 


of the transverse momentum phase space are shown in Fig. for pairs from data and from 
mixed events. The rapidity region covered by the selected pairs essentially extends forward of 
y > 0.5. One observes that only for “Si”- 1 -Si the SSFMs rise faster for the data than for the 
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Figure 13: The second factorial moment F 2 {M) in transverse momentum space for: (a) C+C (window 
of analysis [285,314] MeV), (b) Si+Si (window of analysis [300.9,304] MeV) and (c) Pb+Pb (window 
of analysis [285, 286] MeV) systems. The full triangles represent the moments of NA49 data while the 
open triangles the moments for the corresponding mixed events. 



Figure 14: The combinatorial background subtracted moments AF 2 in transverse momentum space 
for: (a) C+C, (b) Si+Si and (c) Pb+Pb systems. The line in (b) shows the result of a power-law fit 
for > 2000 with exponent 0.33. 


mixed events. The combinatorial backgronnd snbtracted moments AF 2 {M) are plotted versus 
in Fig. M, At larger values of the results are consistent with being constant for “C”+C 
and Pb+Pb whereas one hnds an increase for the “Si”+Si system. Here a power law function 
provides a good £t (y^/dof ~ 0.3) with an exponent $2 = 0.33 ± 0.04 where the error was 
estimated by exploiting the subsample method. The extracted exponent indicates a signihcant 
intermittency effect, but is smaller than the expectation for the CP of ^2 = 0.67. This might 
well be a consequence of the difficulty of isolating the tt+tt" pairs from the a decays. 

Recently another analysis method was developed to estimate critical exponents of fluctua¬ 
tions arising from the existence of a CP. This technique studies hnite size scaling of the particle 
source size parameters as obtained from Bose-Einstein interferometry analysis ED]. The results 
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Figure 15: Scaled variance u of the multiplicity distribution of charged particles versus fj-B for the 
1 % most central Pb+Pb collisions and inelastic p+p reactions for 0 < y < Ubeam (assuming the 
pion mass). Full symbols show results of NA49 m, open symbols preliminary measurements of 
NA61/SHINE [761 [77j. 

based on RHIC and LHC measurements were interpreted as a possible indication of the CP, 
however located at a value of /i^ beyond the range accessible in the SPS energy range. 
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2. Fluctuations of charged particle multiplicity 

The signature of the CP is expected to be primarily an increase of multiplicity fluctua¬ 
tions na which are usually quantihed by the scaled variance u = ((iV^) — {N)‘^)/{N) of the 
distribution of particle multiplicities N produced in the collisions. The measure u is intensive, 
i.e. it is independent of the system volume in statistical models within the GCE formulation. 
However, u is sensitive to the unavoidable volume fluctuations |19]. Therefore the measure¬ 
ments were restricted to the 1 % most central collisions. This selection is based on the energy 
deposited by beam spectator nucleons in the forward calorimeter. Although this tightly con¬ 
strains the number of projectile participants, small fluctuations of the number of the target 
participants remain (see model calculations of Ref. EH)- 


Results for ca of charged particles in Pb-|-Pb collisions (NA49 [75]) are shown in Fig. 15 versus 
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Figure 16: Scaled variance u of the multiplicity distribution of charged particles versus for the 
1 % most central Pb+Pb collisions and inelastic p+p reactions for 1.0 < y < ybeam (assnming the 
pion mass). Full symbols show results of NA49 m, open symbols preliminary measurements of 
NA61/SHINE [761177j . Curves illustrate the effect of a critical point m- 


and compared to preliminary NA61 resnlts from p+p interactions [THIITT]. Evidently, the 
energy dependence of u for the forward hemisphere (which represents almost the fnll acceptance 
of NA49 and NA61) is smooth withont significant maxima. Note that at high SPS energies u 
measnred in inelastic p+p interactions is signihcantly larger than the one for central Pb+Pb 
collisions. The interpretation of this effect is nnder discnssion [771 [78]. 


Figure 16 presents the scaled variance u of charged particles in the rapidity interval 1.0 < 
y < ybeam wliere the azimuthal acceptance of the detector is high and uniform at all energies. 
Again no signihcant irregularities are observed. 

Fluctuations induced by the CP are expected to be stronger at low transverse momenta 
Pt m Therefore uj was also calculated for pt < 0.3 GeV/c and < 0.5 GeV/c. The results are 


plotted in Fig. 17 |75| and turned out to be very similar. 

The effects otf a hypothetical critical point at pb = 360 MeV are illustrated by the curves 
Assuming .^ = 3 fm the value of u was estimated to increase by 0.5 (respectively 


in Fig. 16 


0.25) for all charged (negatively or positively) charged particles [I5l[79| with respect to the value 


^ ys was obtained from statistical model fits to yields of different particle types at the various collision ener¬ 
gies [Uj. It is a monotonically decreasing function of collision energy 
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Figure 17: Scaled variance oj of the multiplicity distribution of negatively charged particles versus 
in the 1 % most central Pb+Pb collisions for 1.0 < y < yteam (assuming the pion mass) [75]. left: 
Pt < 0.3 GeV/c. right: pT < 0.5 GeV/c. Full symbols show data, open symbols are results from the 
UrQMD model. 


expected for the background fluctuations. The limited acceptance of the detector is expected to 
reduce the increase by a factor of the order of 0.6 eg. Guided by the considerations of Ref. [38] 
on the region over which the effects of the CP increase the fluctuations, a parameterization by a 
Gaussian function in ps was chosen with ~ 30 MeV. Evidently, the data do not support 


a maximum as might be expected for a CP (see solid curves in Fig. 16). 

NA49 also obtained u for smaller size nuclei at the top SPS energy of 158R GeV [81]. The 
results together with those for inelastic p+p collisions from NA49 and NA61 are plotted in 


Fig. Interestingly, there may be an indication of a maximum for medium size nuclei. 

A new identiheation procedure (the identity method [H21 IBS]) was developed which allows 
to determine the second and third moments of the multiplicity distribution when the particle 
identiheation is not unique but can only be done on a statistical basis. Applying this method 
NA49 and NA61/SHINE determined the scaled variance of the multiplicity distribution of 
identihed protons, kaons and pions in inelastic p+p and 3.5% most central Pb+Pb collisions. 


The results are shown in Fig. 1+ As in the case of oj for unidentihed charged particles no 
indication of the CP is found. It was pointed out that higher moments of the multiplicity 
distributions are more sensitive to effects of the CP [IH]- The STAR collaboration at RHIC 
performed such a study for the net-proton multiplicity in central Au+Au collisions at energies 
in the range ^/sNN = 7.7 - 200 GeV [El], but also found no evidence for the CP. 
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Figure 18: Scaled variance oj of the multiplicity distribution of charged particles versus the number 
of wounded nucleons N\y in inelastic p+p (1.1 < y < 2.6) and the 1 % most central “C”+C, “Si”+Si 
and Pb+Pb collisions at 158A GeV (1.0 < y < yteam)- Full symbols show results of NA49 [75], open 
symbols NA61/SHINE [76l[77]. 

3. Fluctuations of the net eharge 


Fluctuations of the net charge were originally studied in collisions of heavy nuclei in an 
attempt to hnd evidence for a deconhned phase. The NA49 collaboration concluded that the 
measurements at CERN SPS energies [85] were not sensitive to the initial fluctuations in the 
QGP since they get masked by the effects of resonance decays. The search for the CP has 
rekindled interest in this observable because the CP might enhance net-charge fluctuations 
when the system freezes out in its vicinity. 

Fluctuations of the net charge were quantihed by the measure $ proposed in Ref. ISB and 
dehned in Eq. [T^ For the case of charge fluctuations x is taken to be the electric charge q and 
the measure is called <Fg. 

The distribution of net-charge Q = Xlili Qi ^^al and mixed events is shown in Fig. 


20 


where 


the sum runs over the N particles of the individual events. One observes a clear narrowing 
effect due to charge conservation, which needs to be corrected. In a scenario in which particles 
are correlated only by global charge conservation (CCC) the value of is given by 




g,GCC 


= ^r^-1, 


(19) 
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Figure 19: Scaled variance uj of the multiplicity distribution of protons, kaons, and pions versus 
nucleon-nucleon cms energy yJs^N for the 3.5 % most central Pb-|-Pb collisions (NA49, full symbols) 
and inelastic p-|-p reactions (NA61, open symbols). Large fluctuations for pions in central Pb+Pb 
collisions at high energies are likely due to volume fluctuations. Results are preliminary m- 



Figure 20: The distribution of the net-charge for central Pb-|-Pb collisions at 158A GeV (solid line) 
and the corresponding distribution obtained for mixed events (dotted line) in the maximum rapidity 
interval Ay = 3. 

where 

p (Vol.) 

{^ch) tot 

with {Nch) and {Nch)tot being the mean charged multiplicity in the detector acceptance and in 
full phase space (excluding spectator nucleons), respectively. In order to remove the sensitivity 



30 





















0.1 


O 


-0 


- 0.1 


- 0.2 


-0.3 


* 20AGeV 
^ 30 AGeV 
■ 40 AGeV 


80 AGeV 
160 AGeV 
GCC 


if 


0.1 

c 

o 

^ 0.05 


I 


Hi 




-0.05 


- 0.1 


0.1 0.2 0.3 0.4 

<Nch>/<Nch>tot 


* 20 AGeV • 80 AGeV 
’ 30 AGeV * 160 AGeV 
■ 40 AGeV 



1 T T , , 

: } . 

; ^ V s “ 

I Y Y Y 


Mil? 


3 

Ay 


Figure 21; Left: The dependence of the measure on the fraction of accepted particles for central 
Pb+Pb collisions at 20^-158T GeV. The dashed line shows the dependence expected for the case 
when the only source of particle correlations is global charge conservation. Right: The dependence of 
on the width of the rapidity interval Ay for central Pb+Pb collisions at 20A-158A GeV. Note 
that experimental points for a given energy are correlated as the data used for a given rapidity interval 
are included in the broader intervals. 


to GCC the measure is defined as the difference: 


A*, = *, - *5,GCC 


( 21 ) 


By construction, the value of A$g is zero if the particles are correlated by global charge conser¬ 
vation only. It is negative in case of an additional correlation between positively and negatively 
charged particles, and it is positive if the positive and negative particles are anti-correlated. 

Figure demonstrates that the observed net-charge fluctuations are close to the expecta¬ 
tion for global charge conservation for all the selected bands of rapidity (the charged particle 
multiplicity in the band was used as abscissa). The values Ad)^ stay small at all SPS energies. 


As shown by the energy dependence of A<I)g for two rapidity intervals in Fig. 22 one does not 


observe a peak structure. The STAR collaboration at RHIC performed a study of the higher 
moments of the net-charge distribution in central Au+Au collisions at energies in the range 
^/s~NN = 7.7 - 200 GeV |86] and so far also found no evidence for the CP. 
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Figure 22: The energy dependence of measured in central Pb+Pb collisions at 20A, 30^, 
40^, 80A and 158A GeV for a narrow rapidity interval Ay = 1.2 (left) and a broad rapidity interval 
Ay = 3 (right). 

4- Transverse momentum - multiplicity fluctuations 


Enhancement is also expected for transverse momentum - multiplicity fluctuations when the 
freeze-out occurs close to the CP [T^. For the case of transverse momentum - multiplicity 


fluctuations the fluctuation measures A and S with the IPM normalisation (see Sec. IID) read: 


A[Pt,N] = 


{N)u\p'i 


■[(iV)a;[Pr] - (Pr)a;[iV]] 


and 


where 


and 


j:[Pt,n] = 


{N)uj[pi 


-[{N)uj[Pt] + (PtHN] - 2({PtN) - {Pt){N))] , 


uj[Pt] — 


n;[M] = 


(^1) - (Pt)^ 

(Pt) 

(Al2) _ (iV)2 

W) 


( 22 ) 


(23) 


(24) 


(25) 


are the scaled variances of the two fluctuating extensive event quantities Pt, the sum of the 
absolute values of transverse momenta pti and Al, the number of particles, respectively. The 
quantity u:\pt\ is the scaled variance of the inclusive pt distribution (summation runs over all 
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particles and all events) 


r 1 Pt- Pt^ 

‘^w = 


(26) 


There is an important difference between A[Pt,N] and Tj[Pt,N]. Only the first two mo¬ 
ments: {Pt), {N), and {P^), are reqnired to calculate A[Pt, A^], whereas S[Pt, N] includes 
the correlation term {Pt ■ N). Thus the measures A[Pt,N] and Tj[Pt,N] can be sensitive to 
specihc fluctuations in different ways. Both measures are dimensionless and have a common 
scale required for a quantitative comparison of fluctuations of different, in general dimensional, 
extensive quantities. The values of A and S are equal to zero in the absence of event-by- 
event fluctuations and equal to one for fluctuations given by the model of independent particle 
production. 


The measure is related to the quantity S (see Eq. 14): 

^PT = Vpt'pjIpt] \/I:[Pt, TV] - 1 


(27) 


Results on the dependence of on {^/sJm) in central Pb-|-Pb (NA49 [HIj) and inelastic 


p-|-p collisions (NA61/SHINE preliminary [76l|77]) are plotted in Fig. 23 (top). The measure¬ 


ments are compared to expectations for the CP (solid curves in Fig. 23 top [80]) which were 
obtained in a similar manner like the predictions for a; under the assumption that the increase 
of at the CP amounts to 10 MeV/c for a correlation length of .^ = 3 fm. However, more 
recent theoretical estimates [51] found much less sensitivity of pt fluctuations to the CP. 

Some results on for charged particles from central Au-|-Pb collisions were published by 
the CERES experiment [90] at the SPS for beam energies of 40A, 80A and 158A GeV. The 
results in the pseudo-rapidity acceptance of the experiment (2.2 < p < 2.7) are 1.1 ± 0.4, 
2.3 ± 0.8 and 3.3 ± 0.7 MeV, respectively for the 5 % most central collisions with systematic 
uncertainties of the order of 1.5 MeV. In order to account for a possible change of mean px at 
different beam energies, CERES dehned a dimensionless measure, the ’’normalised dynamical 
fluctuation” S. 


■^Pt- 


(Tt 


, , V pT,dynl 

= sgn{ap^^^y^) ■ — 


-‘Pt — prp^dyn) —— • 

Figure [M| shows that there is no signihcant energy dependence of this measure. 

Measurements by NA49 for different size nuclei at the top SPS energy of 158A GeV are 


shown in Fig. 23 (bottom). As found for multiplicity fluctuations there may also be a maximum 
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Figure 23: Measure of transverse momentum - multiplicity fluctuations of charged particles. Top: 
versus fiB for the 7.2 % most central Pb+Pb collisions (full symbols, NA49 |87j) and inelastic p+p 
reactions (open symbols, NA61/SHINE [76ll77| i. Bottom: versus the number of wounded nucleons Ny/ 
in central C+C, Si+Si and Pb+Pb collisions at 158A GeV (NA49 [88]) and inelastic p+p interactions 
(NA61/SHINE preliminary |76l I77j l. Results are for cms rapidity 1.1 < y < 2.6 assuming the pion 
mass. Curves illustrate the effect of the critical point |80j . 


of transverse momentum fluctuations in medium-size nuclei. 


The corresponding results on Ti[Pt,N] [T^ are presented in Fig. 25 As expected from the 
close relation with they indeed show behaviour consistent with that of Results on 


A[Pr, N] [T^ are shown in Fig. 26 At present there are no predictions for the effect of the CP 
in these observables. 

Finally NA61/SHINE results |7Sl [ZZ] from the two-dimensional scan in system size and 
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Figure 24; The fluctuation measure T,p^ as function of ^ysNN and of /is |90]. The full circles show 
CERES results (after short range correlation (SRC) removal) in central Au+Pb collisions at 40, 80, 
and 158A GeV/c in the pseudo-rapidity range 2.2 < rj < 2.7. The brackets indicate the systematic 
errors. Also shown is the STAR result [89] at -s/snn = 130 GeV which is not corrected for SRC. 
Resnlts and statistical errors from rqmd and URQMD calculations (with re-scattering) are indicated 
as solid and dashed lines, respectively. 

collision energy are presented in Fig. The data come from inelastic p-|-p interactions and 
centrality selected Be-I-Be collisions. No indication of of the ’’critical hill” is observed. 

Data on Ar-|-Sc collisions at 13A, 19A, 30A, 40A, 75A and 150A GeV/c are already recorded 
by NA61/SHINE. They may lead to the discovery of the critical point of strongly interacting 
matter as possibly snggested by the hrst indications seen by the NA49 experiment and discussed 
in this review. 
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Figure 25: Measure Ti\Pt, -/V] (= of transverse momentum - multiplicity fluctuations of charged 

particles. Top: versus fiB for the 7.2 % most central Pb+Pb collisions (full symbols) and inelastic 
p+p reactions (open symbols). Bottom: versus the number of wounded nucleons Nw in inelastic 
p+p and central ”C”+C, ”Si”+Si and Pb+Pb collisions at 158A GeV. Results are for cms rapidity 
1.1 < y < 2.6 assuming the pion mass. (NA49 and NA61/SHINE [76l 177] !. 
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Figure 26: Measure A[Pt, A^] ( = of transverse momentum - multiplicity fluctuations of charged 

particles. Top: versus for the 7.2 % most central Pb+Pb collisions (full symbols) and inelastic p+p 
reactions (open symbols). Bottom: versus the number of wounded nucleons Nw in inelastic p+p and 
central C+C, Si+Si and Pb+Pb collisions at 158A GeV. Results are for cms rapidity 1.1 < y < 2.6 
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Figure 27: Preliminary NA61/SHINE results on Ti[Pt,N] and A[Pt,N] for Be+Be collisions of 
several collision centralities and energies compared with the corresponding results from the energy 
scan with inelastic p+p interactions [761 [TT]. 
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IV. CONCLUSIONS AND OUTLOOK 


The continuing search in nucleus-nucleus reactions for the maximum of fluctuations predicted 
for the critical point of strongly interacting matter has not yet turned up hrm evidence in 
collisions of heavy nuclei either at the CERN SPS (central Pb-|-Pb collisions) or in the RHIC 
BES program (Au-|-Au collisions). 

It is intriguing that both the fluctuations of quantities integrated over the full experimental 
acceptance (event multiplicity and transverse momentum) as well as the bin size dependence of 
the second factorial moment of pion and proton multiplicities in medium-sized Si-|-Si collisions 
at 158A GeV/c possibly suggest critical behaviour of the created matter. 

These results provide strong motivation for the ongoing systematic scan of the phase diagram 
by the NA61/SHINE experiment at the SPS (see Fig. and the continuing search at the 
Brookhaven Relativistic Hadron Collider. 
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